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LOW TEMPERATURE COEFFICIENT OF RESISTANCE AND 
HIGH GAGE FACTOR IN BERYLLIUM-DOPED SILICON 

By J. B. Robertson and M. A. Littlejohn* 

Langley Research Center 

SUMMARY 

The gage factor and resistivity of p-type silicon doped with beryllium was studied 
as a function of temperature, crystal orientation, and beryllium doping concentration. It 
was shown that the temperature coefficient of resistance can be varied and reduced to 
zero near room temperature by varying the beryllium doping level. Similarly, the mag- 
nitude of the piezoresistance gage factor for beryllium -doped silicon is slightly larger 
than for silicon doped with a shallow acceptor impurity such as boron, whereas the tem- 
perature coefficient of piezo resistance is about the same for material containing these 
two dopants. These results are discussed in terms of a model for the piezoresistance o 
compensated p-type silicon. 


INTRODUCTION 

Semiconductor strain gages are nearly 100 times more sensitive to strain than 
metal foil gages and have excellent response times. These semiconductor gages, how- 
ever, have the disadvantage of having large temperature coefficients of resistivity and 
temperature coefficients of gage factor. This problem has been countered primarily by 
use of compensating circuits and by adjusting the impurity concentration in the material 
for minimum temperature effect. Also, electron irradiation of silicon gages (refs. 1,2, 
and 3) has been used to substantially reduce the temperature coefficient of resistivity and 
in some cases to reduce the temperature coefficient of gage factor in these gages. 

Based on the results for electron-irradiated silicon (refs. 1,2, and 3) and the 
physics of the carrier compensation process, one would expect that the thermal charac- 
teristics of silicon strain gages could be varied by any impurity which has a deep energy 
level in silicon. (See ref. 3.) The primary requirements are that the Fermi energy at 
the desired temperature must lie near the energy level in question and that the impurity 
must be present in a large enough concentration to cause the majority carrier concentra- 
tion to vary with temperature near the temperature of interest. Once these requirements 

North Carolina State University. 


2 



are satisfied, then the thermal variations of resistivity and piezoresistance will be deter- 
mined primarily by the impurities present in the material, the condition of charge 
neutrality, and the effects of strain on the energy levels and the semiconductor band 
parameters. 

It has been shown that the temperature coefficient of resistivity can be easily 
reduced to a zero value near room temperature by diffusing beryllium into high purity 
p-type silicon. (See ref. 4.) This indicates that the majority carrier concentration does 
vary with temperature in this temperature range. Also, it is known that beryllium has 
two acceptorlike centers in silicon whose ground state energies are at 0.145 eV and 
0.190 eV above the valence band edge. (See ref. 5.) These results were deduced from 
infrared absorption spectroscopy. However, when the temperature dependence of the Hall 
coefficient is studied, these two levels appear to be one acceptor level at approximately 
0.168 eV above the valence band edge (ref. 6). Thus, for typical values of mobility in 
p-type silicon (ref. 7), the Fermi energy would be located near this "average level" for 
resistivities of approximately 1 ohm -cm. Resistivities near this value are easily 
achieved by beryllium diffusion. Thus, it appears that beryllium satisfies the necessary 
requirements for varying the thermal properties of p-type silicon strain gages. 

This paper presents the results of piezoresistance and resistivity measurements 
on beryllium -doped silicon gages which show that these gages can be made which have a 
zero temperature coefficient of resistivity at room temperature, a gage factor over 200, 
and a temperature coefficient of gage factor equal to that of boron -doped silicon gages, 
which are the most frequently used type of silicon gages. 

SYMBOLS AND ABBREVIATIONS 
Eg e beryllium activation energy 

Ej Fermi energy 

f Fermi probability function 

G-p gage factor 

gage factor due to band changes only 
k Boltzmann constant, 8.616 X 10~ 5 eV/K 

N a acceptor concentration 
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N Be 

N V 

P 

Pi 

Po 

q 

R o 

R(e) 

T 

A Be 

AE V,eff 

AE Vi 

AE f 

A f> A V,eff 

e 

M 


beryllium concentration 

valence band effective density of states 

hole concentration 

hole concentration in ith valence band 
hole concentration in unstrained material 
electronic charge 
resistance for zero strain 
resistance as a function of strain 
temperature, kelvins 

change in beryllium energy level due to strain 

effective shift in valence band edge 

shift in individual valence bands 

change in Fermi energy with strain 

deformation potentials (see eqs. (15)) 

strain 

carrier mobility 


p i 

P'eff P 

j u mobility of holes in silicon in ith band 

a hole mobility 


v 


Poisson’s ratio 



p electrical resistivity 

o conductivity, (ohm-era)"^ 

TCGF temperature coefficient of gage factor 

TCR temperature coefficient of resistivity 

EXPERIMENTAL TECHNIQUES 

The material used was obtained by thermal diffusion of beryllium into high purity, 
single -crystal, p-type silicon (p > 100 ohm-cm). Silicon wafers, 2.5 cm in diameter and 
0,10 cm thick, were plated on both sides with beryllium using vapor deposition. The 
plated wafer was sandwiched between two clean silicon wafers and placed in a furnace in 
a helium environment. The samples were heated at 1300° C for 1 hr. After diffusion, 
the resistivity was approximately 0.2 ohm-cm p-type, which represents the solubility 
limit of beryllium in silicon. To obtain resistivities greater than 0.2 ohm-cm the beryl- 
lium concentration was reduced by heating the silicon at 1100° C until the desired resis- 
tivity was reached. Resistivities as high as 7 ohm-cm could be easily achieved. 

The wafer was mechanically lapped to a thickness of about 0.025 cm. The resis- 
tivity profile over the surface was examined by using a four -point probe. The variation 
in the resistivity over the region of the wafer from which samples were cut was typically 
20 percent. After final lapping the crystal orientation was determined by X-ray diffrac- 
tion and by optical goniometer. 

Bridge -shaped strain-gage samples were cut from the wafers with an ultrasonic 
cutter. The samples were then chemically polished and ohmic contacts formed with 
evaporated aluminum. The aluminum was alloyed at 600° C for 10 minutes in an argon 
atmosphere. Immediately after alloying, the aluminum contacts were nickel plated to 
enhance the uniformity of the contact. The final active dimensions of the strain -gage 
samples were approximately 0.51 cm by 0.191 cm by 0.018 cm. 

The strain gages were mounted to stainless -steel cantilever beams having dimen- 
sions 15.2 cm by 2.54 cm by 0.254 cm. First, the area of the beam to which the gage was 
to be mounted was sandblasted and chemically cleaned. An initial coat of strain-gage 
epoxy was applied to the beam and cured for 1 hr at 250° C. Next, an additional coat of 
epoxy was applied, the sample was then fixed to the beam, and again the epoxy was cured 
at 250° C for 1 hr. 
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The cantilever beams were mounted in a cryostat, the bottom section of which is 
illustrated in figure 1. The temperature of the samples could be varied from about 90 K 
to 400 K, although most measurements were taken from 220 K to 37 5 K. Strain was 
applied to the sample by deflecting the beam with a micrometer which was thermally 
insulated from the beam. Both tensile and compressive strains could be applied. 

The calculation of the strain in the beryllium -doped gage is very important in deter- 
mining the actual magnitude of the piezo re si stance gage factor. (This problem has been 
discussed in reference 8 for an end- loaded cantilever beam such as that used in the pres- 
ent investigation.) The approach taken herein was to first examine several (four in all) 
commercial gages, both metal and semiconductor, to determine the linearity of the sys- 
tem in figure 1. The results were quite good in that for each gage tested the normalized 
change in resistance as a function of the normalized deflection yielded a straight line with 
a slope of unity to within less than 5 percent. The same results were obtained for strains 
up to about 10' 3 cm/cm (calculated from the manufacturers’ specified gage factors). 

Also, the strain levels deduced from the manufacturers’ specified gage factors were in 
agreement with each other to within about 5 percent. However, these strains were con- 
sistently lower than the theoretical strain calculated from the beam dimensions and its 
deflection, as much as 15 percent lower. Since one of the purposes of the present study 
was to compare piezoresistance gage factors of beryllium -doped gages with those of 
boron-doped gages, the strain used in the calculations was the theoretical strain in an 
end- loaded cantilever beam. It is felt that this will provide a reasonably lower limit on 
the experimentally determined piezo re si stance gage factors. 

After a sample was mounted in the cryostat and contact wires were attached, the 
cryostat was assembled and evacuated to a pressure of 1.33 mPa (1CT^ mm Hg). All 
measurements were taken under vacuum in the dark. The resistance was measured by 
passing a constant current from a current calibrator through the sample and measuring 
the voltage across the sample with a calibrated digital voltmeter. The accuracy of the 
resistance measurement was better than 0.1 percent. The temperature of the samples 
was controlled with a temperature controller and measured with a copper-constantan 
thermocouple. The accuracy was at least 1 percent, and caution was taken to insure 
that thermal equilibrium was established before any measurement was taken. 

The nature of the epoxy bond was checked by examining the resistance change as a 
function deflection in both tension and compression up to strains of about 10 " 3 cm /cm. 
Also, the gage was cycled between a strained and unstrained state many times and finally 
was left in a strained state for at least an hour before any thermal cycles were begun. 

Any indication of an appreciable built-in strain or creep in the bond resulted in the sam- 
ple being discarded. The room -temperature gage factor was then measured at a strain 
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of about 4.5 X 10'^ cm /cm. Next, the piezo re si stance was studied as a function of tem- 
perature. After the sample was returned to room temperature, the gage factor was 
measured again to make sure that the room temperature value had not changed due to 
degradation of the epoxy bond. 


EXPERIMENTAL RESULTS 


Transport Properties 

The electrical transport properties of a material are of considerable interest in 
considering it for strain-gage applications. Thus, since one objective of this work was to 
compare strain gages of beryllium -doped silicon to boron-doped silicon, the transport 
properties of beryllium -doped sillicon were studied. 

The electrical resistivity p of any material is given by 


P 


1 _ 1 

cr nq/i 


( 1 ) 


where a is the conductivity, q is the electronic charge, n is the carrier concentra- 
tion, and /i is the carrier mobility. For a p-type semiconductor, 


where p is the hole concentration and ju p is the hole mobility. 

The dependence of resistivity on temperature will be determined by the dependence 
of the carrier concentration and the carrier mobility on temperature. 


P(T) “ qp(T) Hp(T) (3) 

Carrier mobility decreases with increasing temperature for the temperature range of 
interest. For holes in lightly doped silicon, the dependence is approximately p p cc T“^*^ 
Carrier concentration in semiconductors increases with increasing temperature. In 
extrinsic (doped) silicon at temperatures less than 250° C, the carrier concentration is 
determined by thermal ionization of impurities. At sufficiently low temperatures, the per 
centage increase in carriers with temperature ^/p is greater than the percentage 
decrease in mobility with temperature ^/u p , and consequently, the TCR will be negative 
As the temperature increases, more carriers are ionized, and the supply of un-ionized 
impurities diminishes; thus, the increase of carriers is slowed. At sufficiently high 
temperatures, the percentage increase in carriers is less than the percentage decrease 
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in mobility, and, consequently, the TCR will be positive. At some intermediate tempera- 
ture, therefore, the TCR must be zero, and the resistivity pass through a minimum. 

This occurs at the temperature where the increase in carriers is equal to the decrease 
in mobility. This condition is obtained mathematically by setting the derivative of 
equation (3) equal to zero, which gives 

dp/dT _ ~ dM p/ dT ( 4 ) 

P 

The ionization energy of boron in silicon is 0.045 eV, and most of the boron is 
ionized at temperatures well below room temperature. As a result, the resistivity min- 
imum of boron-doped silicon occurs between -175° C and -75° C, The effective ionization 
energy of beryllium in silicon is ~0.168 eV. Therefore, a higher temperature is required 
to thermally ionize a fraction of beryllium impurities than is required to ionize the same 
fraction of boron impurities. As a result, the resistivity minimum in beryllium -doped 
silicon occurs at a higher temperature. 

Figure 2 shows a plot of Hall mobility as a function of temperature for samples 
having hole concentrations at room temperature ranging from 2.55 x 10 15 cm -3 to 
1.65 x 10 17 cm" 3 . Also shown for comparison is the curve for T’ 2 ’ 7 dependence on 
temperature which is found for lightly doped p-type silicon (ref. 7). Figure 3 illustrates 
the dependence of the hole mobility on the hole concentration. The solid curve is based 
on calculations from an empirical relation given in reference 9, which was derived from 
experimental data reported in references 10, 11, and 12. Thus, the mobility values of 
beryllium -doped silicon and their temperature and carrier concentration dependence are 
in good agreement with reported values for silicon doped with shallow impurities such as 
boron. 

Figure 4 shows the experimentally determined resistivity as a function of tempera- 
ture for several different beryllium doping levels. The room -temperature resistivity 
and hole concentration are indicated for each curve. It is noted that the temperature at 
which zero TCR (i.e., the minimum in each curve) occurs can be selected from a wide 
temperature range by adjusting the beryllium doping level. For a room -temperature 
resistivity of about 1.0 ohm-cm, the TCR is approximately zero. At room temperature, 
1.0-ohm -cm boron-doped silicon has a TCR of s=0.75%/K. Shown in figure 4 is a T 2 - 7 
temperature dependence which would correspond to the variation of the resistivity of 
boron-doped silicon for doping levels below about 5 x 10 16 cm" 3 , and for temperatures 
above about 0° C. (See ref. 7.) 
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Piezoresistance 


The experimental results presented herein are in terms of the longitudinal piezo- 
resistance gage factor which is given by 



R(e) - R 0 


V 


( 5 ) 


where R(e) is the resistance under strain, R Q is the resistance for zero strain, and 
e is the strain along the gage axis. This is the customary sensitivity factor for strain 
gages and is exact if the resistance is linear with strain. The linearity of resistance 
with strain was examined for each sample tested. Essentially no deviations from linearity 
were observed for strains up to 10 " 3 cm/cm. All gage factor measurements were taken 
at a strain of between 4.5 x 10 cm/cm to 4.7 x 10" 4 cm/cm. Small changes in the 
strain level resulted from slight changes in the location of individual samples on the 
cantilever beam. 

Figure 5 shows measured values of the piezoresistance gage factor taken at room 
temperature as a function of the sample resistivity for various crystal orientations. 

The solid lines and squares are reported data for boron-doped silicon. The solid lines 
(ref. 13) represent the highest gage factor reported in the literature for boron-doped 
silicon. 

In general, the gage factor of beryllium -doped silicon is larger than for boron -doped 
silicon. Again, it must be pointed out that the data points represent conservative values, 
since the strain used in equation (5) is the theoretical strain at the surface of the deflected 
cantilever beam. By using the strain values obtained from a commercially available 
strain gage, the data points in figure 5 would be at least 10 percent higher. 

The temperature dependence of the gage factor was measured for several gages. 

The data for five of the gages are shown in figure 6 and are representative of the temper- 
ature dependence of all the gages studied. The change in gage factor is approximately 
30 percent per 100° C and the TCGF is only very slightly dependent on the doping level; 
this is also true for boron -doped silicon. (See ref. 1.) The solid curve in figure 6 
represents the temperature dependence of gage factor in 1.0-ohm- cm boron -doped silicon 
and, for all practical purposes, is the same as that for the beryllium -doped gages. 


DISCUSSION 


Two factors regarding the piezoresistance need to be discussed. First, the data 
in figure 5 indicate that the difference in the room -temperature piezoresistance gage 
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factor between beryllium -doped and boron-doped silicon is nearly independent of the 
doping level. Second, the data of figure 6 indicate that the temperature coefficient of 
piezoresistance gage factor for beryllium -doped silicon is almost the same as that of 
boron-doped silicon. 

For small strains, the piezoresistance gage factor can be approximated as (ref. 2) 


G„ = l + 2,-I^ 

F CT de 


( 6 ) 


where v is Poisson's ratio and a is the conductivity, which in the strained state is 


= q £ M iPi 


( 7 ) 


where ^ is the mobility and is the concentration of the holes in the ith valence 
band. The conductivity can be written as 


cr=q(i ef£ P 


( 8 ) 


where 


eff 




0) 


and 


p = 2 p i 


Thus, 


°F = 1 + 2 * - JT 


1 d/J eff 1 dp 


eff 


de P de 


( 10 ) 


The term 


1 d ^eff 


is, to first order in the strain, the influence on the gage factor due 


M e ff de 

to band edge changes with strain alone. Thus, 


G - G^ . — 
U F “ U F p de 


(ID 


Here Gj?^ is the gage factor due to band changes and would be the gage factor 


of 
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boron-doped silicon in the temperature range where all the boron centers are ionized and, 
thus, cannot change the hole concentration with strain. 

Now, the hole concentration under strain can be written 


P - P Q exp 



( 12 ) 


where AE^ is the change in the Fermi energy with strain, p 0 is the unstrained hole 
concentration, and AEy is an effective shift in the valence band edge with strain 
given by 



Here AEyj is the shift of the individual band edges with strain. Equation (12) assumes 
no changes in the effective masses of the individual bands with strain. (See ref. 14.) 

At small strains, the shifts in the energy levels are taken to be of first order in the 
strain, (ref. 14) so that for small strains 


1 dp = A f " A V,eff 
P de kT 


(14) 


where A^ and Ay are "deformation potentials," which satisfy 
AE f = 

AE V,eff = A V,eff e 

The hole concentration can also be written as 



(15) 


p = fN 


N 


Be 


Be 


( E Be 

+ ex H^tr 


- Ej 


( 16 ) 


where f is the Fermi probability function, Ng e is the beryllium concentration, Eg e 

is the beryllium activation energy, and E^ is the Fermi energy. Differentiating equa- 
tion (16) with respect to e and using equations (15) gives another expression for 
as follows: 
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( 18 ) 


Combining equations (17) and (14) to eliminate Aj gives 

1 dp 1 - f A Be ~ A V,eff 

Pdf 2 - f kT 

Thus, the piezoresistance gage factor of beryllium -doped silicon is 


G 


F 



1 - f A Be “ A V,eff 

2 - f kT 


(19) 


An estimate of the change in the beryllium energy level due to strain Ag e 
can be made from the data presented so far. At 400° C for the lightly doped samples 
p _ N + Ng e « Ng e . Thus for the 4.34-ohm-cm samples, Ng e x4x 10*5 cm"3, and 
at room temperature f * 0.8. Also Gp - G^ ~ 20 in the [ill] direction. Thus, 

A Be “ A V,eff _ 12Q 
kT 

Using data given by Wortman et al. (ref. 14) gives Ay e g = -3.3 eV/unit strain. Then 
A Be * -0*2 eV/unit strain. It should be pointed out that this calculation assumes that 
the beryllium energy level does not ’’split" with strain, but only shifts in energy. 

CONCLUDING REMARKS 

Beryllium -doped silicon strain gages were made which had a zero temperature 
coefficient of resistance at room temperature. The temperature at which zero tempera- 
ture coefficient of resistance occurs can be adjusted over a considerable range above 
and below room temperature by adjusting the beryllium concentration. The beryllium- 
doped gages have slightly higher gage factors, lower temperature coefficients of resist- 
ance and comparable temperature coefficients of gage factor than the standard boron- 
doped silicon gages. The results imply that similar properties may be obtained by 
doping with other "deep” impurities such as indium whose ionization energy in silicon 
is 0.16 eV. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., Febuary 21, 1974. 
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Figure 1.- Schematic of bottom section of 
cryostat used to take piezoresistance 
and resistivity measurements. 
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Figure 2.- Hall mobility as a function of temperature 
for beryllium -doped silicon. 
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Figure 3.- Hole mobility as a function of hole 
concentration for beryllium -doped silicon. 
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Figure 4.- Resistivity as a function of temperature 
for beryllium -doped silicon mounted on beam. 
Numbers given with each curve represent: 
room-temperature resistivity; room- 
temperature carrier concentration. 
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Figure 5.- Room -temperature piezoresistance gage 
factor as a function of resistivity for beryllium - 
doped silicon and boron-doped silicon. 



Figure 6.- Normalized piezoresistance gage factor 
as a function of temperature for beryllium - 
doped silicon. 
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